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Abstract 

The astrophysical S-factor for proton-proton fusion, Sn(E), is obtained with the nuclear matrix 
element analytically calculated in pionless effective field theory. To the third order, the zero-energy 
result Sn(0) and the first energy derivative S u (0) are found to be (3.99 ±0.14) x 10~ 25 MeV b 
and Sn(0) (11.3 ±0.1) MeV -1 , respectively; both consistent with the current adopted values. 
The second energy derivative is also calculated for the first time, and the result S 11 (0) = 
Sn(0) (170 ± 2) MeV -2 contributes at the level of 0.5% to the fusion rate at the solar center, 
which is smaller than 1% as previously estimated. 
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I. INTRODUCTION 



As the trigger of stellar hydrogen burning, the proton-proton (pp) fusion reaction, 
p+p —> d+e + +V e , plays a fundamental role in astrophysics. Its reaction rate, conventionally 
encoded by the astrophysical S-factor Sn(E), with E the energy in the center of mass 
frame, is thus an important input in studies such as stellar evolution and solar neutrinos. 
For such a pivotal reaction, the value of Su(E) at the center of the Sun (with temperature 
~ 1.55 x 10 7 K which yields E ~ a few keV), unfortunately can not be determined or reliably 
extrapolated by terrestrial experiments. This is because they have to be performed at much 
higher energies to overcome the Coulomb barrier and get sensible statistics. Therefore, one 
has to rely on theory for predictions. 

Ever since the first proposal of the pp chains by Bethe and Critchfield the nuclear 
transition amplitude of the pp fusion, together with other solar fusion cross sections have 
been extensively studied. Those results were extensively reviewed first in Ref. j^J (SFI) and 
then most recently in Ref. js| (SFII). For solar fusion, the temperature is much lower than 
typical energy scales in nuclear physics, thus only the first few terms in the expansion of 
Sn(-E) around E = is needed for solar models. Currently, the recommended value of 
Sn(0) in SFII has an error of 1%, which results in a 1% error in the pp fusion rate. The 
recommended value of S 11 (0) also has an error of 1%, which only results in a less than 
0.1% error in the fusion rate. S 11 (0) was not given in SFII. However, Bahcall and May {4] 
estimated its contribution to the rate to be ~ 1%, comparable to the overall error in Su. 
Thus, SFII recommended a modern calculation of S n (0) be undertaken. 

In this work, we calculate Sn(.E') using the nuclear effective field theory (EFT) with pions 
integrated out. One of the major ingredients, the square of the orbital matrix element A 2 (E) 
which Sn(-E') is proportional to, is extracted from the cross section of an analogue process: 
v e + d — > p + p + e~, analytically computed in Ref. Q]. 

This pionless EFT is applicable for low energy processes with the characteristic momen- 
tum p much smaller than the pion mass m n 6|-[8[ , which is the case for solar pp fusion. In this 
theory, pions are integrated out. All the nucleon-nucleon interactions and two-body currents 
are described by point-like contact interactions with a systematic expansion in powers of 
p/rUn- A close analogy of this theory is the Fermi theory of four fermion contact interactions. 
The one- and two-body contributions both depend on the momentum cut-off but the sum 
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does not. In pionless EFT, there is only one two-body current (with coupling in every 
weak interaction deuteron breakup process to next-to-next-to-leading (NNLO) in the p/m n 



expansion {5]. This two-body current is a Gamow- Teller operator. The other two-body 
currents are either missing due to vector current conservation or the matrix elements are 
suppressed because of the orthogonality of the initial and final state wave functions in the 
zero recoil limit. This means the universal number L^a encodes the two-body contributions 
for all low energy weak deuteron breakup processes, and it takes just one measurement to 
calibrate all the processes. This feature is also seen in the other complimentary approaches 



to the solar fusion such as potential models [9|, hybrid version of EFT 10J, and pionless 
EFT with dibaryon (see SFII for more details). 

II. 8u(E) FROM PIONLESS EFT 

The astrophysical S-factor, S(E), for a nuclear reaction at kinetic energy E is related 
to the reaction cross section by 

a (£) = ^e- 2 ^. (1) 

The rapid- varying energy dependence of cr(E) due to the Coulomb barrier is mostly ac- 
counted for in the exponential term which is controlled by the Sommerfeld parameter 

where a = 1/137.036 is the fine structure constant and m p = 938.272 MeV is the proton 
mass. 1 

Using the same convention and inputs as SFII, we write the S-factor for the pp fusion as 
Sn(-E) = 6vr m p a In 2 3— (— ) 777^ , (3) 

where 7 = yj 2 fi np E d = 45.70 MeV (/i np = 469.459 MeV is the proton-neutron reduced 
mass and E^ = 2.224573 MeV is the deuteron binding energy) is the deuteron binding 
momentum; gy = 1 and qa = 1.2695 are the usual Fermi and axial- vector coupling constants, 
respectively; and (ft) +_^ + = 3071 sec is the ft value for superallowed Fermi transitions. 



1 Wc use the units h = c = 1. 
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The energy dependence of Sn(E) is determined by two terms: (i) A(E), the orbital matrix 
element, which is proportional to the nuclear transition matrix element, and (ii) f^(E), the 
phase space factor in this nuclear (3 + process. We note the conventional way of separating the 
radiative correction in Sn(E): the long-distance (so-called "outer") part, which is process- 
dependent, is included in the phase factor f£L(E) (annotated by a superscript "if), while the 
short distance (so-called "inner") part, which is process-independent, is taken into account 



by the product of (gA/gv) 2 an d l/(/*)o+->o+- For a detailed account, see Ref. [11]. 



A. Orbital Matrix Element A(E) 

At very low energy, the pp fusion predominantly goes from the 1 5'o partial wave state to 
the deuteron state ( 3 S\ with some 3 Lh mixture) through the spatial axial current operator. 



In pionless EFT, it takes the form 12] 



A -=9A N ^ T - akN 

+L 1>A [(A T P fc A)t(jV T P~AO + h.c] +■•• , (4) 

where N denotes the nucleon field; r~ = (ti — ir 2 ); V = V — V; Pk = t 2 o 2 ak/y/8; 
P = t 2 t~ o"2/v / 8; all a's and r's are the Pauli matrices for spin and isospin, respectively. 
The coupling constant L\ tA of the leading axial two-body current appears at next-to-leading 
order (NLO) and there is no new two-body current contributing until next-to-next-to-next- 



to-leading order (N LO). The orbital matrix element A(p), with p = ^m p E, is then related 
to the nuclear matrix elements of A k by 

\(d;j\A;\pp}\=g A C v ip)^^A(p)8i, (5) 
where j denotes the deuteron \d) polarization and 

is the square of the Sommerfeld factor. The effect of the deuteron recoil, with momentum 
q < 0.4 MeV, is suppressed by a factor q 2 /'-/ 2 < 10 -4 and hence can be neglected. In the zero 
recoil limit the vector matrix element between pp and d vanishes because the wave functions 
are orthogonal. 
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InRef. [5|, various (anti)neutrino deuteron breakup processes were computed in pionless 
EFT up to N 2 LO with analytic expressions. The needed hadronic matrix element for the 
pp fusion can be extracted from the result of the similar process v e + d — > p + p + e~ . The 
square of the orbital matrix element can be simplified and expressed as 

A 2 (P) = 5-^5 t^- [Fiip) + Kip) + G' 2 {p)\ + 0(e 3 ) , (7a) 

where pd = 1.764 fm is the effective range parameter for deuteron, and the corresponding 
energy-dependent structure functions are 



2m P 7p 

7T (p 2 + 7 2 ) 



CHp)e 



4r,(p) tan-l(£) 



2 V 



Fjp) =-Im 

7T 



B 2 (p) (A^(p) + e A (p) + e 2 Ai(p)) 



(7b) 
(7c) 



e — 1m 

7T 



m 2 7 




2 ix z 



2 -^B' (p)L 1A (A^(p) + eA (p)) 



L'IaA.^p) 



(7d) 



The complex integral 



£o(p) = m p 



(2vr) 3 fc 2 + 7 2 



p 2 — k 2 + id 



(8) 



<5->0 



The components of the pp scattering amplitude in the 1 S'o channel, Ai(p), are found 
to be 



A-i (p) 
4>(p) 
>4i(p) 



47T 
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[1/a + a m p 


H(V)} ' 
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[1/a + a m p 


^(r/)] 2 


7T 
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[1/a + a m p 


H( V )]* 



(9) 



where a = —7.82 fm and ro = 2.79 fm are the scattering length and effective range parameter 
for 1 Sq, respectively; and the function 



d{irf) 2ir] 



ln(i 77) 



(10) 
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depends on momentum via the Sommerfeld parameter r](p). The renormalization scale /i- 
independent coupling constants L[ A is defined as 



{p~l) 
m p Co-i l 



Li, 



m p 

2 IT 



pd 



(11) 



(P ~7) S 

and when being evaluated at the pion-mass scale, \x = m n « 140 MeV, the low energy 
constants Cq-i = — 3.77fm 2 and 62,-2 = 7.50 fm 4 . 2 

Note that we keep the small expansion parameter e explicit in Eqs. ( |7aT - l7d[) to make the 
order of each term transparent. For evaluation, one has to make a series expansion in e up 
to the maximum order at which the result is valid, and then set e = 1 in the end. 

At zero energy, our result of A 2 (0) can be written in a compact form 

1 



A 2 (0) 



where the function 



[e* - 7 a(l-xe x E 1 ( X )) ~ ei 2 aL lA f + 0(e 



(12) 



Ei{x) 



00 e -t 

dt — 
t 



(13) 



yields a value of 1.4655 with \ — & m p/l — 0.1498. Up to the order of e 2 , this result is in 
agreement with Ref. 12| in which A(0) is worked out to the order of e 4 ; and the constant 
Li : a defined therein is exactly the same as our L x A . This result is also consistent with 
Ref. 13] by ignoring the L' lA term, and Ref. jl^] by a redefinition of L' lA . 



B. Phase Space Factor fJ*(E) 



The phase space factor f(E) in a nuclear /3 process is conventionally written as (see, e.g. 



Ref. 



I* 



1 rQ+E 

f(E) = — dW F(±Z, W)p e W(Q + E- Wf , (14) 

where m e = 0.511 MeV is the electron (and positron) mass; W and p e = \JW 2 — ml are 
the relativistic energy and momentum of the emitted electron or positron; Z is the charge 
of the final nuclear state with the sign "±" being assigned to /3 T emission, respectively; and 



As we take the q 2 /j 2 — > limit in F\. F^, and Bo, and re-define G2 and Li,a by some normalization 
factors, we add a prime to remind readers about these changes from Ref. 
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Q is the difference of the total rest mass of the initial and final nuclear states. The Fermi 
function 

(15) 

where 70 = (1 — Z 2 a 2 ) 1 / 2 and v e 



[r(2 7o + i)] 2 ' 

^1 — Z z a z ) L/z and v e = ±Z a W/p e take both the relativity of electron and the 
finite size of nucleus (through a spherical radius R) into account. 

For pp fusion, Z = 1, R = 2.1402 fm, Q = 0.420 MeV, and the outer radiative correction 



evaluated to be 5°^ 



1.62% 
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so the "corrected" phase space factor becomes 



L R P (E) 



[l + SZ^fppiE) 



pp 



0.144 1 + 9.04 



MeV J 



30.7 



The linear term in E is in excellent agreement with [4|. A 0.1% error is assign to /^(0) in 
SFII. An error of the same size is assigned to f^JE) in this work. 



VMeV J J 



+ 0(E 3 ). (16) 



III. RESULTS AND DISCUSSION 



To obtain a numerical result of A. 2 (p), the last piece of information we need is the value 



of Li t A- Using the reactor u e d breakup data, Li.a is determined to be (3.6 ±4.6) fm 3 18l. Il9|. 
and using the data of solar neutrino deuteron breakups through charged and neutral currents, 
and v e e~ elastic scattering, L\ a is determined to be (4.0 ±6.3) fm 3 191 ] . Both fits have taken 
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17j |. Treating these two fits in 



into account the radiative corrections computed in Refs. 
two nucleon systems as independent and using the same average scheme as in Ref. [18| . we 
obtain 

L M = (3.4±3.7)fm 3 , (17) 

where a — 0.3fm 3 shift of the central value is introduced due to updating qa from 1.267 to 
1.2695. This range of Li a is consistent with the nave dimensional analysis value \Li a 



6fm d [20J. One can further constrain this two-body current by the tritium beta decay 
which is a three nucleon system. This is carried out in potential model 9J and hybrid 
EFT [lo|, and both approaches yielded the determination of Sn(0) with 1% accuracy. 
Expanding A 2 in powers of E, 

MeV J \MeV J 



A 2 (£)=A 2 (0) 1 + d 



0(E 3 



(18) 



we find, by using the central value of L^a, the following expansions in e for A 2 (0), c\ and 
c 2 : 



A 2 (0) oc 1 + 0.057e + 0.054e 2 + 0.022e 3 + ■ ■ • , 
ci oc 1 - 0.112e + 0.049e 2 + 0.020e 3 + ■ • ■ , 

c 2 oc 1 - O.OOOe + 0.057e 2 + 0.023e 3 + ■ • • . (19) 

Although these results are only good up to 0(e 2 ), the above 0(e 3 ) terms indicate ~ 2 — 3% 
corrections. This is what one typically expects for the pionless EFT, which has a small 
expansion parameter ~ 7/771^ ~ 1/3. Therefore, we assign 3% higher order corrections to 
A 2 (0), Ci and c 2 . 

Up to 0(e 2 ), the full results are 

A 2 (0) = 7.01 (1 ± 3% ± 3%) , 
ci = 2.24 (1 ± 0.7% ± 3%) , 

c 2 = 34.2(1 ±0.1% ±3%), (20) 

where the first errors are from the uncertainty in and the second errors are from higher 
order (0(e 3 )) corrections. 

The two sources of errors are correlated in the sense that when two near-threshold weak 
deuteron breakup processes have the same kinematics, the two processes will be governed 
by the same Gamow- Teller matrix element such that higher order effects can be largely 
absorbed by shifting the value of L^a- Indeed, the two nucleon processes used to constrain 
L\^a and the pp fusion have similar kinematics - even though not exactly the same, so the 
combined error on A 2 (0) should be between 3% (100% correlated) to 4.2% (not correlated, 
added in quadrature). We will assign the combined error to be 3.6% which is slightly larger 
than the 3% error adopted in SFII with pionless EFT determination. Note that in SFII, a 
0.9% error is assign to A 2 (0) by using tritium beta decay rate to constrain the two body 
current. Here we just use the constraint of Eq. ffTTl) available in the two nucleon sector and 
study its implication in Su(E). 

It is remarkable that c\ and c 2 are quite insensitive to L\^. This is because for near 
threshold pp fusion, the energy dependence is dominated by the initial state pp scattering 
which alone gives C\ = 3.09 and c 2 = 36.2 without L^a dependence. c\ and c 2 then receive 
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—28% and —6% Li^-dependent corrections, respectively, to reach the values of Eq. ([2~0]) . 
This explains why the errors due to the uncertainty of L^a in c\ and are 23% and 4% of 
that in A 2 (0). This result along with the combined error assignment discussed above leads 
to error reduction in S n (0)/Sn(0) and S 11 (0)/Sn(0). 

Combining our results of A 2 (E) and f pp (E), we found the S-factor Sn(.E') for the low- 
energy pp fusion yields 

Sn(0) = (3.99 ±0.14) x l(T 25 MeVb, 
S' n (0) = Sn(0) (11.3 ± 0.1) MeV" 1 , 

Si'^0) = Sn(0) (170 ± 2) MeV" 2 . (21) 

Our central value of Sn(0) is the same as the pionless EFT value in SFII while the error is 
slightly larger as discussed above. As for S n (0), our result is consistent with the adopted 
value in SFII: s'® m (0) = (11.2 ± 0.1) Sn(0)MeV _1 , which was obtained by Bahcall and 



May [4|. Finally, we also report a value of S 11 (0), which was not computed before, with ~ 1% 
accuracy. At the center of the sun, this higher-order term has a contribution at the level of 
0.5% to the fusion rate, which is smaller than ~ 1% that was previously estimated [3j, |4j. 
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